fsamp = 5000;
fcuts = [800 1000];
mags = [1 0];
devs = [0.04 0.01];
[n,Wn,beta,ftype] = kaiserord(fcuts,mags,devs,fsamp);
hh = fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');
freqz(hh)
